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ABSTRACT 

We present archive ESO VLT photometric and spectroscopic data of the Ultraluminous X-ray source 
NGC 1313 X-2. The superb quality of the VLT images reveals that two distinct objects, with R magni- 
tudes 23.7 and 23.6, are visible inside the Chandra error box. The two objects, separated by 0.75", were 
unresolved in our previous ESO 3.6 m+EFOSC image. We show that both are stars in NGC 1313, the 
first a BO-09 main sequence star of ~ 20M©, while the second a G supergiant of ~ IOMq. Irrespectively 
of which of the two objects the actual counterpart is, this implies that NGC 1313 X-2 is a high mass 
X-ray binary with a very massive donor. 

Subject headings: galaxies: individual (NGC 1313) — stars: individual (NGC 1313 X-2) — X-rays: 
binaries — X-rays: galaxies 



1. INTRODUCTION 

First revealed by Einstein, the population of ultralu- 
minous X-ray Sources (ULXs) has increasingly grown up 
in the last decade mainly thanks to the observations of 
ROSAT (e.g. Colbert & Ptak 2002), XMM-Newton (e.g. 
Foschini et al. 2002a) and Chandra (e.g. Liu & Breg- 
man 2005a; Swartz et al. 2004). About 150 ULXs are 
included in the recent Chandra catalogue of Swartz et al. 
(2004). These point-like sources have X-ray luminosities 
Lx ~£i 10 39 erg s" 1 , in excess of that of a ~ 10M Q com- 
pact object accreting at the Eddington limit. Variability 
in the X-ray flux on timescales of months is observed in 
about half of the i?OS^4T ULXs with multiple observations 
(Colbert & Ptak 2002), while about 5-15% of the Chandra 
ULXs show variability during a single observation (average 
exposure time ~ 40 ks, Swartz et al. 2004). 

For several sources with sufficiently good statistics, the 
best fit to the X-ray spectrum is obtained with a two- 
component model, a soft multicolor disk (MCD) black- 
body plus a power law. Some ULXs show typical temper- 
atures of the MCD component 5-10 times lower than those 
of Galactic X-ray binaries. The high luminosity, the very 
soft thermal component (if it represents the emission from 
a cool accretion disk) and the variability suggest that these 
sources may be powered by accretion onto an Intermediate 
Mass Black Hole (IMBH) of 100-1000 M Q . Nevertheless, 
many of the ULX properties can be explained if they do 
not emit isotropically (King et al. 2001) or are dominated 
by emission from a relativistic jet (e.g. Kaaret et al. 2003). 
In this case, they may harbor stellar mass BHs and may 
be similar to Galactic black hole binaries. 

Multiwavelength observations are definitely a powerful 
tool to investigate the nature of ULXs. Radio emission, 
when present, gives important clues on the geometry, en- 
ergetics and lifetime of ULXs (Kaaret et al. 2003; Miller et 
al. 2005). Optical follow-ups are crucial to identify ULX 
counterparts and clean up the population from the signif- 
icant contamination of background AGNs and interacting 



SNe (Foschini et al. 2002b; Masetti et al. 2003; Swartz et 
al. 2004) . Up to now only a very small number of ULXs 
have been convincingly associated with stellar objects of 
known spectral type (e.g. Liu et al 2002, 2004; Kaaret et al. 
2004). All these ULXs are hosted in star- forming regions 
and their optical counterparts have properties consistent 
with those of early type O-B stars. Some ULXs are also 
associated with extended optical emission nebulae (Pakull 
& Mirioni 2002). 

NGC 1313 X-2 is a prototypical ULX (see Miller et al. 
2003; Zampieri et al. 2004; Turolla et al. 2005, and ref- 
erences therein). With a luminosity Lx ~ 10 40 erg s _1 in 
the 0.2-10.0 keV band, it is a good candidate for harboring 
an IMBH (M > 100M Q ). Such an option is corroborated 
by the presence of a very soft X-ray spectral component 
(T ~ 200 eV) which points to a compact object of mass 
definitely larger than those of Galactic Black Hole candi- 
dates. Moreover, the object exhibits X-ray variability on 
a timescale of months. 

On the basis of a 19 ks Chandra exposure and accu- 
rate astrometry of field objects, Zampieri et al. (2004) 
(Z04 hereafter) derived the source position with an un- 
certainty of 0.7" (RA=03:18:22.34, DEC=-66:36:03.7; la 
confidence level) . Inside the Chandra error box a faint op- 
tical candidate was found on a R band image taken with 
the ESO 3.6 m telescope in January 2002. 

Here we present a follow-up study of the optical counter- 
part of NGC 1313 X-2, based on photometric archive data 
obtained with the ESO VLT telescope and a reconsidera- 
tion of our previous photometry of the ESO 3.6 m+EFOSC 
image. 

2. OPTICAL OBSERVATIONS: DATA REDUCTION AND 
ANALYSIS 

We analyzed archive ESO VLT+FORS1 images (BVR) 
and spectra of NGC 1313 X-2 taken between December 
2003 and January 2004 (Program ID 072.D0614). The 
observations are listed in Table 1. Data reduction was 



1 Department of Astronomy, University of Padua, Vicolo dell'Osservatorio 2, 1-35122 Padova, Italy; mucciarelli@pd.astro.it 

2 INAF- Astronomical Observatory of Padua, Vicolo dell'Osservatorio 5, 1-35122 Padova, Italy; zampieri@pd.astro.it, falomo@pd.astro.it 

3 Department of Physics, University of Padua, Via Marzolo 8, 1-35131 Padova, Italy; turolla@pd.infn.it 

4 Department of Physics and Mathematics, University of Insubria, Via Valleggio 11, 1-22100 Como, Italy; treves@mib.infn.it 



2 



Mucciarelli et al. 



performed using standard IRAF 5 procedures. The images 
were astrometrically calibrated using 29 GSC2 stars. The 
calibration uncertainty, tested with GSC2 stars not used 
for the calibration, is ~ 0.3". Finally, for each band the 
images were combined and cleaned from cosmic rays. Fig- 
ure 1 shows the combined R and B images. On the same 
night a short exposure of the standard PG 0231+051 Lan- 
dolt field (Landolt 1992) was also taken in each band. 
Aperture photometry (1" radius) was performed on the 
combined images. The instrumental magnitudes were then 
calibrated with the Landolt standard stars in the Bessel- 
Cousins system (see Patat 2003, for extinction coefficients 
and color terms). 

Z04 give the R magnitude of a number of objects around 
NGC 1313 X-2 and of the proposed counterpart (object C 
in their paper). The latter was close to the limit of de- 
tcctability on their image and appeared as a single object. 
Thanks to the higher resolution of the VLT image, in the 
R and V exposures we are able to resolve object C in 
two distinct point sources, CI and C2. Both are inside 
the Chandra error box (see Figure 1). Object C2 is not 
detected in the B band frame. Magnitudes, colors and as- 
trometric positions of the two candidate counterparts, CI 
and C2, and of objects A, B, and D (following Z04) are 
reported in Table 2. The photometric errors are the 2a 
statistical errors on the measurements with the different 
Landolt standards. For object C2, we quote an upper limit 
to the B band magnitude using the plate limit (B = 25.2). 

The total magnitude of object C1+C2 is R = 23.2±0.05. 
We note that this measurement is not in agreement with 
the estimated magnitude of object C obtained in our pre- 
vious ESO 3.6m image (R = 21.6; see Z04). In fact, re- 
analyzing the old image we found that the stellar back- 
ground used for subtraction was underestimated, leading 
to a measurement ~ 0.3 mag brigther. We also found a 
mistake in the adopted exposure time of the image that 
causes an overestimate of another ~ 0.9 mag. The revised 
magnitude of object C in the ESO 3.6m image is now found 
to be R = 22.9 ±0.2, in reasonable agreement (within 2cr) 
with the VLT measurement. 

In addition to the images, we also analyzed four 
VLT+FORS1 spectra (A c =5900 A, A/AA=440 at A c ) of 
objects C1+C2 taken in different nights (Table 1). The 
slit (1") was oriented to include object D. After perform- 
ing standard reduction, all spectra were sky subtracted, 
wavelcnght calibrated through comparison lamp exposures 
and flux calibrated using standard star spectra obtained 
in the same night. In these VLT+FORS1 spectra the two 
sources (CI and C2) are not spatially resolved. The 2D 
spectrum taken on 15 January 2004 is shown in Figure 2. 
Nebular emission lines of [Oil] A 3727 A, H 7 , FL3, [OIII] AA 
4959-5007 A, [OI] A 6300 and 6364 A, H Q , [Nil] A 6583 A 
and [SII] AA 6717-6731 A are clearly detected. Note that 
this is the first detection of a [Oil] line from this nebula. A 
one dimensional spectrum was extracted over an aperture 
of 2.2" centered on object C1+C2 from each of the four 
combined spectra. Two nebular spectra were extracted 
from different 1" apertures, eastward and westward of the 
source position and adjacent to it. The two spectra were 
then averaged and the resulting spectrum subtracted from 

5 IRAF is distributed by the National Optical Astronomy Observato 
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that of object C1+C2. All these spectra, taken on Jan- 
uary 15th, are shown in Figure 3. The nebula-subtracted 
source spectrum shows no evident emission or absorption 
lines. Residuals are present in coincidence with some neb- 
ular lines (expecially [OIII] and H Q ), with an upper limit 
to the equivalent width of ~ 30 A. In particular the resid- 
ual flux in the [OIII] line is a non negligible fraction of the 
nebular flux. This appears to be caused by an increased 
emission of the nebular line around the position of object 
C1+C2. It is not clear if this is simply induced by a change 
in the rather irregular spatial profile of the nebular line or 
by a variation of the physical conditions produced by the 
presence of the nearby ULX. Finally, marginal evidence 
of an excess in emission may be seen at 4686 A, corre- 
sponding to Hell emission, but the line is not statistically 
significant. 

3. RESULTS 

The superb quality of the VLT images reveals that two 
distinct objects, CI and C2, are visible inside the Chan- 
dra error box of NGC 1313 X-2 in the R and V bands. 
From the astrometric positions reported in Table 2, we 
infer a separation of 0.75" and a position angle (C2 with 
respect to CI) of ^131° . Note that the two objects were 
unresolved in our previous 3.6 m+EFOSC image. 

The possibility that either CI or C2 may be a back- 
ground AGN appears very unlikely. In fact, no statistically 
significant emission line at wavelengths longer than H Q is 
observed in the optical spectrum nor any other feature 
that may be identified with a highly redshifted emission 
line. Furthermore, from the R magnitude of objects CI 
and C2, and the X-ray flux of the last XMM-Newton ob- 
servation (fx ~ 2 x 10~ 12 erg cm~ 2 s _1 ; Z04), we obtain 
a X-ray to optical flux ratio fx/ fa ^ 2000. This is more 
than an order of magnitude larger than the typical value 
for an AGN (< 100). The possibility that object C2 is 
a highly obscured AGN is not very convincing cither. In 
fact, to produce the same X-ray flux but a fx/ fa ratio 
of ~ 100, the R magnitude should be ~ 20. Therefore 
the additional reddening required to make object C2 look 
like an obscured AGN is Ar ~ 3.5 mag, corresponding to 
a column density ~ 9 x 10 21 cm~ 2 (Bohlin et al. 1978), 
much larger than that inferred from X-ray spectral fits 
[~ (3±g;|) x 10 21 cm" 2 ]. Unless absorption in NGC 1313 
is characterized by an extinction law quite different with 
respect to that in our Galaxy, we are led to rule out an 
identification with an (obscured) AGN. 

A photometric analysis of a large sample of field stars 
(> 30) was performed on the three BVR images. A color- 
color plot of these objects is shown in Figure 4. Almost all 
of them have stellar colors, apart from two on the lower 
right part of the diagram that are bluer than ordinary 
stars. Given their apparent magnitude, objects earlier 
than spectral type F (including CI) are stars in NGC 
1313, while the others may be Galatic foreground stars 
or belong to NGC 1313. However, object C2 cannot be a 
Galatic foreground star because its absolute visual magni- 
tude would be My <~ 9, too large to be consistent with its 
colors. 

Within the photometric errors, the colors of object CI 

;, which are operated by the Association of Universities for Research in 
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appear to be consistent with those of a A3-09 I or a A2- 
BO V star, while those of C2 with a G8-G7 I star (see e.g. 
Cox 2000; Braddley 1982). Unfortunately, the optical con- 
tinuum does not provide useful information for assessing 
the spectral type because the light from both objects con- 
tributes to it. Observationally, the slope of the continuum 
can be caracterized by a power law, A -18 . The absence 
or extreme weakness of the Hell A 4686 A emission line in 
the optical spectrum suggests that X-ray irradiation is not 
dominant. Taking Galactic absoption into account and as- 
suming Ay ~ 0.3 (Cardelli ct al. 1989 exctinction law with 
Rv — Av/Eb-v = 3.1 has been used throughout), the 
unreddened colors of object CI are (V — R)o — —0.2 ± 0.2 
and (B — V)o = —0.2 ± 0.2, consistent with those of a 
B8-0 I or AO-05 V star (see Figure 4). For object C2 it is 
(V - R) = 0.4 ± 0.2 and (B - V) > 1.0, consistent only 
with a G4 I star (see again Figure 4). Recently, Liu et al 
(2005b) performed a 6.4 m Magcllan/Baade observation of 
the field around NGC 1313 X-2 and found a I = 23.3 mag 
object in coincidence with the position of C1+C2 (that 
appear unresolved in their I frame). Assuming that the 
flux in the / band originates mainly from the redder object 
C2, we then obtain (R — I) = 0.3 ± 0.2, consistent with 
our tentative spectral classification. 

At the distance of NGC 1313 (d = 3.7 Mpc; Tully 1988) 
the V magnitudes of CI and C2 (reported in Table 2) 
translate into the absolute magnitudes My ~ —4.6 and 
~ —4.1, respectively. Comparing these values with the ab- 
solute magnitudes of main sequence and supergiant stars 
(e.g. Cox 2000; Braddley 1982), we find that the observed 
value is consistent only with a BO-09 main sequence star 
for CI, while it is consistent with a G4 supergiant of type 
lb for C2. Therefore, we conclude that both CI and C2 are 
stars in NGC 1313, with CI an early type main sequence 
star and C2 a supergiant. The bolometric luminosities of 
the two objects are ~ 3 x 10 38 erg s _1 and ~ 2 x 10 37 erg 
s _1 , respectively. 

Given the density of objects in the field of view (~ 
50 — 100 arcmin -2 ), a significant fraction of which are 
supergiants in NGC 1313, the probability that CI or C2 
fall by chance inside the 2a Chandra error box is not neg- 
ligible (~ 0.1). However, the chance occurrence of two 
objects, separated by only 0.7", inside the X-ray error box 
is ~ 5 x 10~ 3 , sufficiently small to be considered rather 
unlikely. Actually, if both CI and C2 are stars in NGC 
1313, a physical association may not be unplausible (the 
distance corresponding to the apparent separation on the 
sky is ~ 10 pc). Therefore, we conclude that the ULX 
is most probably physically associated to either object CI 
or C2. Irrespectively of which of the two objects is the 



actual counterpart, NGC 1313 X-2 appears to be a high 
mass X-ray binary with a very massive donor star. 

We now briefly consider what are the consequences if ei- 
ther CI or C2 are the counterparts of the X-ray source. A 
BO-09 main sequence star has an initial mass of ~ 2OM0. 
In this respect our analysis essentially confirms the original 
suggestion by Z04, who proposed that the optical counter- 
part of NGC 1313 X-2 may be an O type main sequence 
star in NGC 1313. If the colors are affected by the binary 
interaction, the estimated mass may vary somewhat. A 
~ 2OM donor star could easily provide the mass trans- 
fer rate required to fuel the accreting black hole through 
Roche-lobe overflow during the main sequence phase if the 
orbital separation is w 1 AU. In these conditions, the mass 
transfer would be stable and the source persistent (Pa- 
truno et al. 2005; Patruno & Zampieri, in preparation). 

If object C2 is the counterpart, the nature of the system 
remains unchanged. In this case X-ray irradiation may be 
significant and give a non-negligible contribution to the 
optical emission. The mass corresponding to a G4 super- 
giant is ~ IOMq. The same caveat discussed for CI about 
the possibility that the colors and mass estimate are af- 
fected by binary interaction applies also in this case. The 
mass transfer rate provided by such a donor star throug 
Roche-lobe overflow is certainly adequate also for large 
orbital separations. Wind accretion may also be a viable 
alternative. 

At the present stage, there seems to be no compelling 
argument to rule out any of the two objects CI and C2 
as possible counterparts of NGC 1313 X-2. In both cases, 
NGC 1313 X-2 is a high-mass X-ray binary with a very 
massive donor. In order to distinguish which of the two 
objects is the counterpart, further observations are needed. 
In particular, we aim at optical spectroscopy under opti- 
mal seeing conditions and with the slit along the line con- 
necting CI and C2, which would allow us to check the 
existence of typical lines of X-ray binaries (like He II A 
4686 A). Theoretical modelling (binary evolution calcula- 
tions) may also help distinguish the counterpart. It will 
also clarify the properties of the accretion flow and the 
compatibility with the IMBH or stellar black hole inter- 
pretations. 

Based on observations made with ESO Telescopes at 
the Paranal Observatories, obtained from the ESO/ST- 
ECF Science Archive Facility. This work has been par- 
tially supported by the Italian Ministry for Education, 
University and Research (MIUR) under grant PRLN-2004- 
023189. We thank Alessandro Patruno for interesting dis- 
cussions on the properties of IMBH binary sistems. 
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Fig. 1.— R (left) and B {right) VLT+FORS1 images of the field around NGC 1313 X-2 (30"x20"). The circle is the 2a Chandra error-box 
(1.4"). In the R frame, the counterpart is clearly resolved in two point sources, CI and C2. 
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Fig. 2. — Two-dimensional spectrum (VLT+FORS1) of objects C1+C2. The slit (1") is oriented to include object D. The wavelength 
intervals are 3600-5200 A (left) and 6200-7100 A (right). 



Table 1 

Log of the VLT+FORS1 observations. 



Obs. type 


Date 


Filtcr/Grism(A c ) 


Exp. time (s) 


Seeing (") 


Image 


2003-12-24 


B 


840x2 


1.0 


Image 


2003-12-25 


V 


600x2 


0.8 


Image 


2003-12-24 


R 


500x2 


0.8 


Spec. 


2003-12-22 


5900 A 


1300x2 




Spec. 


2003-12-24 


5900 A 


1300x2 




Spec. 


2003-12-30 


5900 A 


1300x2 




Spec. 


2004-01-15 


5900 A 


1300x2 





Table 2 

astrometric positions, magnitudes and colors of the sources around ngc 1313 x-2 (see figure 1). 



Source 


RA 


DEC 


B 


V 


R 


B-V 


V-R 


A 


03:18:21. 97±0.05 


-66:36:06.4±0.3 


23.5±0.15 


21.7±0.05 


20.6±0.05 


1.8±0.15 


l.liO.l 


B 


03:18:21.57±0.05 


-66:36:00.8±0.3 


22.4±0.15 


22.7±0.05 


22.5±0.05 


-0.3±0.15 


0.2±0.1 


CI 


03:18:22.26±0.05 


-66:36:03.3±0.3 


23.5±0.15 


23.6±0.15 


23.7±0.15 


-0.1±0.2 


-0.1±0.2 


C2 


03:18:22.36±0.05 


-66:36:03.8±0.3 


>25.2 


24.1±0.15 


23.6±0.15 


>1-1 


0.5±0.2 


D 


03:18:20. 96±0.05 


-66:36:03.6±0.3 


20.3±0.15 


18.9±0.05 


18.1±0.05 


1.4±0.15 


0.8±0.1 
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Fig. 3. — From top to bottom: VLT+FORS1 spectrum of 15 January 2004 (3700-8900 A), extracted from an aperture of 2.2" centered on 
the unresolved object C1+C2; average spectrum of the nebula, extracted from two 1" apertures located eastward and westward of the source 
position and adjacent to it; the nebula-subtracted spectrum of object C1+C2. 
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(V-R). 



Fig. 4. — (B — V)o vs. (V — R)o for a sample of field objects around NGC 1313 X-2, including A, B, CI, C2 and D. Error bars are shown 
at the upper left corner. Measurements were corrected for Galactic exctinction (Ay = 0.3, Eg_y = 0.1, corresponding to a column density 
Nfj = 6 X 10 20 cm -2 ). The arrow indicates the reddening vector corresponding to Eg_y = 0.1. The solid and dashed lines represent the 
colors of supergiant (SG) and main sequence (MS) stars. 



